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Abstract. 

We present three dimensional hydrodynamical simulations aimed at studying the dynamical 
and chemical evolution of the interstellar medium (ISM) in isolated dwarf spheroidal galaxies 
(dSphs). This evolution is driven by the explosion of Type II and Type la supernovae, whose 
different contribution on both the dynamics and chemical enrichment is taken into account. 
Radiative losses are effective in radiating away the huge amount of energy released by SNe 
explosions, and the dSph is able to retain most of the gas allowing a long period (^2 — 3 Gyr) 
of star formation, as usually observed in this kind of galaxies. We are able to reproduce the stellar 
metallicity distribution function (MDF) as well as the peculiar chemical properties of strongly 
O-depleted stars observed in several dSphs. The model also naturally predicts two different 
stellar populations, with an anti-correlation between [Fe/H] and velocity dispersion, similarly to 
what observed in the Sculptor and Fornax dSphs. These results derive from the inhomogeneous 
pollution of the SNe la, a distinctive characteristic of our model. We also applied the model to 
the peculiar globular cluster (GC) uj Cen in the hypothesis that it is the remnant of a formerly 
larger stellar system, possibly a dSph. 

Keywords, hydrodynamics, methods: numerical, stars: abundances, ISM: evolution, ISM: abun- 
dances, galaxies: dwarf, (Galaxy:) globular clusters: individual (a) Cen) 



1. Introduction 

Due to their proximity, the galaxies of the Local Group (see Mateo 1995 and Geisler 
et al. 2007 for a review) offer an unique opportunity to study in details their structural, 
dynamical and chemical properties and to test different theories of galaxy formation. 

Owing to their low metallicity and lack of neutral hydrogen, it was initially believed 
that dSphs are relatively simple objects whose ISM is completely removed by SN II 
explosions after a very short intense star formation period (e.g. Dekel & Silk 1986). 
Doubts about this picture come from the high resolution spectroscopy of several dSphs 
showing a wide range in metallicity. For istance, Shetrone et al. (2001) have observed 
stars in Draco and Ursa Minor with values of [Fe/H] in the range —3 ^[Fc/H]^ —1.5. 
The same authors also found that their observed dSphs have [a/Fe] abundances that 
are 0.2 dex lower than those of Galactic halo field stars at the same metallicity. This 
suggests that the bulk of the stars in these systems formed in gas self polluted by SNe II 
as well as SNe la and that the star formation (SF) must continue over a relatively long 
timescale in order to allow a sufficient production of iron by SNe la (Ikuta & Arimoto 
2002, Lanfranchi & Matteucci 2004, Marcolini et al. 2006; but see also Recchi et al. 2007, 
Salvadori et al. 2008). 
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Figure 1. Left: distribution of the logarithm of the density in three orthogonal planes at ~ 400 
Myr after the beginning of the simulation. Note the inner bubble carved by the SNe explosions 
and the shocks propagating outward. Right: [Fe/H] distribution of the models resembling the 
Draco dSph (left panels) and the ui Cen globular cluster (right panels) together with the cor- 
responding [a/Fe]-[Fe/H] diagrams; the blue line and blue dots represent stars with [a/Fe]5C0.2 
(i.e. affected by SNe la inhomogenous pollution, see text). The lower panels represents the 
assumed SFH for the two models. 

A complex star formation history (SFH) is further suggested by several facts: i) isolated 
low mass dSphs such as Phoenix (Young et al. 2007) and Leo T (de Jong et al. 2008) 
were able to form stars up to 100 Myr ago; ii) the SFHs of dwarf galaxies are strongly 
dependent on their local environment, the fraction of passively evolving galaxies dropping 
from ^70% in dense environments, to zero in the rarefied field (Haines et al. 2007 ); Hi) 
dwarf ellipticals and dSphs cluster around the dominant spirals galaxies, while gas rich 
star forming dwarf Irregulars are found at larger distances (van den Bergh 1994). These 
points highlight the role of the environment (tidal interaction/ram pressure stripping), 
and disfavors a scenario in which the evolution is due uniquely to internal processes. 

Here we briefly report some results obtained by Marcolini et al. (2006, 2007, 2008) 
with a model which turns to be consistent with many properties of the Draco dwarf and, 
with minimal assumptions, to the chemical properties of the peculiar system oj Centauri 
(Marcolini et al. 2007), which is believed to be the remnant of an ancient dSph. 



2. dSphs Model 

Let us consider the reference model of Marcolini et al. (2006) which, although employed 
to study the general characteristics of dSphs, is taylored to explore the evolution of 
the Draco dSph. The simulation starts with the ISM in hydrostatic equilibrium in the 
extended dark matter halo potential well. The amount of initial gas corresponds to the 
cosmological baryonic fraction of the dark matter halo (Mjsm = 0.18Mdm)- The SFH 
is given a priori, assuming that stars form in a sequence of 50 instantaneous bursts 
separated in time by 60 Myr (see Fig. I). We also assume that SNe II explode at a 
constant rate for 30 Myr after the occurrence of each starburst, while the SNe la rate 
follows the prescription of Matteucci & Recchi (2001). Each SN explosion is stochastically 
placed into the galaxy according to its radial probability P{r) — M+(r)/Mi< t tot, where 
M+(r) and M+^ot are the nowadays radial stellar mass profile and total stellar mass, 
respectively (cf. Marcolini et al. 2006 for more details). 

Although the total energy released by the SNe II explosions is larger than the binding 
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Figure 2. Abundance ratio [a/Fe] versus [Fe/H] for 500 sampled stars of the reference model 
at t = 2 Gyr, compared to a dataset of 28 dSph stars (Draco: triangles ; Fornax: crosses; Leo I: 
asterisks; Sextans: diamonds; Sculptor: circles; Ursa Minor: squares) collected from the literature 
(see Marcolini et al 2008 for details). 

energy of the ISM, efficient radiative losses enable the galaxy to retain most of its gas, 
which thus remains available for the aforementioned prolonged SFH. The burst of SNe 
II associated with each stellar burst pushes the bulk of the ISM to the outskirt of the 
galaxy (see Fig. 1, left panel). Once the explosions cease (~30 Myr after each star burst 
episode), the ISM flows back towards the center of the galaxy; when the next burst 
occurs, the gas is pushed outwards again. This oscillatory behavior leads to a rather 
efficient and homogeneous pollution of the ISM by the SNe II ejecta. We note that as 
the SFH has been fixed a priori, it has no direct relation to the gas reservoir within the 
galaxy. However, stars do form during the quiescent periods between bursts, when the 
gas has "settled" , and an a posteriori consistency for the SFH is recovered. For example 
a similar periodic ISM behaviour was recovered in similar low-mass galaxies simulations 
performed by Stinson et al. (2007). We point out that, while the observed dSphs are 
deprived of gas, the galaxy in the present model can not expel the bulk of its ISM by 
internal mechanisms. Following the discussion given in the Introduction, an external 
cause, such as the interaction with the Milky Way (Mayer et al. 2006), must be invoked 
at some evolutive stage of the galaxy to get rid of its ISM. 

Given their lower rate, SNe la do not significantly affect the general hydrodynamical 
behavior of the ISM, but their role is relevant for the chemical evolution of the stars. 
Because of their longer evolutionary timescales, SN la progenitors created in previous 
starbursts continue to explode during the subsequent quiescent periods, when the gas is 
flowing back into the central region. During these periods the higher ambient gas density 
(together with the lower SNe la explosion rate) cause the SNe la remnants to be isolated 
from one another, forming chemically inhomogeneous regions (we refer to these regions 
as "SNe la pockets"). These pockets are "washed out" by successive phases of expansion 
and collapse of the ISM, due to the effects of SNe II, but new pockets form during the 
quiescent phases between consecutive starbursts. At odd with the ejecta of SNe II, SN la 
debris are rich in iron and deficient in a-elements. Thus, stars forming in the SN la 
pockets possess lower [a/Fe] and higher [Fe/H] ratios than those formed elsewhere. 

2.1. Comparison with Local dSphs 

As stressed above, the reference model discussed in Marcolini et al. (2006) was tailored 
to fit the Draco galaxy. In Fig. 1 we show the MDF of this model at the end of the 
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Figure 3. Age-Z and age-[Fe/H] distributions of 1000 sampled stars for the w Cen model. 
Note the metallicity spread among coeval stars. 

simulation (3 Gyr). The maximum value of the distribution occurring at [Fe/H]~-1.5 
and the mean value of ([Fe/H])=-1.65 with a spread of ~1.5 dex are compatible with 
observations (Bcllazzini et al. 2002) . The agreement is larger at the high metallicity tail of 
the distribution which, in our model, is shaped by the stars formed in the SN la pockets. 

As discussed above, the SNe la are also responsible of the [cn/Fe] spread present at 
larger values of [Fe/H] in the observed dSphs, as shown in Fig. 2 (lower panel). From 
this figure it is also apparent that the most metal poor stars have [a/Fe]~ 0.5, typical 
of pure SN II enrichment, while a decrement of ~ 0.2 dex in the plateau is achieved at 
higher metallicities. This drop occurs in our simulations after t ~ 2.0 — 3.0 Gyr to allow 
the cumulative effect of SNe la to be appreciable (top panel of Fig. 2). We thus conclude 
that, at least in our models, only a prolonged star formation history (> 2.0 Gyr) can 
account for the chemical differences between the Galactic halo and the dSphs as shown 
by Shetrone et al. (2001). Such a long SFH for dSphs is also found in the chemical models 
by Fenner et al. (2006) who are unable to reproduce the Ba/Y ratio unless stars formed 
over an interval long enough for the low-mass stars to pollute the ISM with s-elements. 

The inhomogeneous pollution by SNe la discussed above is particularly important 
in the central galactic region, where the SN la rate is greater, and the density of the 
ambient gas is higher. This naturally accounts for a radial segregation of Fe-rich stars in 
the central regions of dSphs. As these stars in our model are also a-depleted, a similar 
radial segregation of a-poor stars should be observed to test our model. We finally stress 
that our model envisages a central depression in the radial distribution of the stellar 
velocity dispersion (Marcolini et al. 2008). This naturally entails two different stellar 
populations with an anti-correlation between [Fe/H] and velocity dispersion, which has 
been observed in the case of the Sculptor dSph (Tolstoy et al. 2004) and the Fornax dSph 
(Battaglia et al. 2006). 

2.2. uj Cen 

The stellar system ui Cen (NGC 5139) is unique among Galactic star clusters in terms 
of its structure, kinematics, and stellar content. It is the only known GC showing a 
clear [Fe/H] spread spanning the metallicity range —1.6 <[Fe/H]< —0.6 (Norris et al. 
1996, Johnson et al. 2008). Recent photometric surveys have revealed the presence of 
multiple sequences in its color-magnitude diagram (CMD), indicating a complex star 
formation history (e.g. Sollima et al. 2005). It is suggested that ui Cen is the nucleus 
of a larger stellar system, possibly a dwarf galaxy, that lost most of its stars and gas 
in the interaction with the Milky Way ~10 Gyr ago (e.g. Bekki & Norris 2006). In this 
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framework Marcolini et al. (2007) focused on the evolution of the central region of their 
dSph model, where the inhomogeneous pollution of the SNe la is particularly effective. 

In Fig. 1 the MDF within the inner 90 pc is plotted assuming a total SFH lasting ~ 1.5 
Gyr as shown in the same Figure (see Marcolini et al. 2007 for more details and about 
the assumptions used for this model). The MDF shows a bimodal structure similar to 
that observed in u> Cen (e.g. Norris et al. 1996), with a maximum at [Fe/H]=— 1.6 and 
a secondary peak at [Fe/H]^ —1.3 accounting for ~ 25% of the cluster's stellar content. 
Comparing the [a/Fe]-[Fe/H] diagram of this model with the one typical of dSphs (c.f. 
Fig. 1), it is possible to note the much larger number of a-depleted stars (see also the blue 
line in the MDF which represents stars with a/Fe^0.2). The diagram is in reasonable 
agreement with the findings of Pancino et al. (2002) and Origlia et al. (2003), who 
find that while the metal-poor and intermediate-metallicity stellar populations of to Cen 
have the expected a-element overabundance observed in halo and GC stars (([a/Fe]) 
~ 0.3-0.4), the most metal-rich population ([Fe/H]~ —0.6) shows a significantly lower 
a-enhancement (([a/Fe]) ~0.1). 

Given the distictive role of the SN la pollution, the [Fe/H] content of this system is not 
simply proportional to the metal content and a large spread in the age-metallicity relation 
is always present (see Fig. 3). These peculiarities have important consequences on the 
cluster CMD. In fact, they reduce by a factor of 40% the large (and still unexplained) 
helium overabundance usually invoked to account for the anomalous position of the blue 
main sequence observed lo Cen. 
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